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Thermodynamic properties of a strongly coupled quark-gluon plasma (QGP) of constituent quasi- 
particles is studied by a color path-integral Monte-Carlo simulations (CPIMC). For our simulations 
we have presented QGP partition function in the form of color path integral with new relativistic 
measure instead of Gaussian one used in Feynman and Wiener path integrals. For integration over 
color variable we have also developed procedure of sampling color variables according to the group 
SU(3) Haar measure. It is shown that this method is able to reproduce the available quantum lattice 
chromodynamics (QCD) data. 



I. INTRODUCTION 



Ph. 

Oji Investigation of properties of the QGP is one of the main ehallenges of strong-interaetion physics, both theoretically 

^ ■ and experimentally. Many features of this matter were experimentally discovered at the Relativistic Heavy Ion 

,^ ' Collider (RHIC) at Brookhaven. The most striking result, obtained from analysis of these experimental data [l|, is 
that the deconfined quark-gluon matter behaves as an almost perfect fluid rather than a perfect gas, as it could be 
expected from the asymptotic freedom. There are various theretical approaches to studying QGP. Each approach 
^ , has its advantages and disadvantages. The most fundamental way to compute properties of the strongly interacting 
matter is provided by the lattice QCD [2|, Q ■ Interpretation of these very complicated Monte Carlo computations 

0^ \ requires application of various QCD motivated, albeit schematic, models simulating various aspects of the full theory. 

^^ . Moreover, such models are needed in cases when the lattice QCD fails, e.g. at large baryon chemical potentials and 

^ ; ' out of equilibrium. While some progress has been achieved in the recent years, we are still far away from having a 

CO . satisfactory understanding of the QGP dynamics. 

^^ ' A semi-classical approximation, based on a point like quasiparticle picture, has been introduced in [^-|8|. It is 

^^ , expected that the main features of non-Abelian plasmas can be understood in simple semi-classical terms without the 
difficulties inherent to a full quantum field theoretical analysis. Independently the same ideas were implemented in 
terms of molecular dynamics (MD) [3 . Recently this MD approach was further developed in a series of works [10, [ill ■ 
The MD allowed one to treat soft processes in the QGP which are not accessible by perturbative means. 

A strongly correlated behavior of the QGP is expected to show up in long-ranged spatial correlations of quarks 
C^ . and gluons which, in fact, may give rise to liquid-like and, possibly, solid-like structures. This expectation is based 
on a very similar behavior observed in electrodynamic plasmas |12l . Il3i . This similarity was exploited to formulate a 
classical nonrelativistic model of a color Coulomb interacting QGP [lO| which was numerically analyzed by classical 
MD simulations. Quantum effects were either neglected or included phenomenologically via a short-range repulsive 
correction to the pair potential. Such a rough model may become a critical issue at high densities, where quantum 
effects strongly affects properties of the QGP. To account for the quantum effects we follow an idea of Kelbg [ij 
rigorously allowing for quantum corrections to the pair potential |27| . To extend the method of quantum effective 
potentials to a stronger coupling, we use the original Kelbg potential in the path integral approach, which effectively 
map the problem onto a high-temperature weakly coupled and weakly degenerate one. This allows one to extend the 
analysis to strong couplings and is, therefore, a relevant choice for the present purpose. 

In this paper we extend previous classical nonrelativistic simulations jlOf based on a color Coulomb interaction to 
the quantum regime. For quantum Monte Carlo simulations of the thermodynamic properties of QGP we have to 
rewrite the partition function of this system in the form of color path integrals with a new relativistic measure instead 
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of Gaussian one used in Feynman and Wiener path integrals. For integration partition function over color variables 
we have developed procedure of sampling color quasiparticle variables according to the group SU(3) Haar measure 
with the quadratic and cubic Casimir conditions. The developed approach self-consistently takes into account the 
Fermi (Bose) statistics of quarks (gluons) . The main goal of this article is to test the developed approach for ability 
to reproduce known lattice data [2| and to predict other properties of the QGP, which are still unavailable for the 
lattice calculations. First results of applications of the path integral approach to study of thermodynamic properties 
of the nonideal QGP have already been briefly reported in [la, [i3- In this paper we have shown that CPIMC is able 
to reproduce the QCD lattice equation of state and related thermodynamic functions and also yields valuable insight 
into the internal structure of the QGP. 

Theoretical treatment and hydrodynamic simulations of the experimentally observed expansion of the fireball con- 
sisting of quarks and gluons at relativistic heavy-ion collisions suggests knowledge not only thermodynamic but also 
transport properties of the QGP. Unfortunately the CPIMC method itself is not able to yield transport properties. 
To simulate quantum QGP transport and thermodynamic properties within unified approach we are going to combine 
path integral and Wigner (in phase space) formulations of quantum mechanics. The canonically averaged quantum 
operator time correlation functions and related kinetic coefficients will be calculated according to the Kubo formu- 
las. In this approach CPIMC is used not only for calculation thermodynamic functions but also to generate initial 
conditions (equilibrium spatial, momentum, spin, flavor and color quasiparticle configurations) for generation of the 
color trajectories being the solutions of related differential dynamic equations. Correlation functions and kinetic co- 
efficients are calculated as averages of Weyl's symbols of dynamic operators along these trajectories. The basic ideas 
of this approach have been published in [l8| . Using this approach we are going to calculate diffusion coefficient and 
viscosity of the strongly coupled QGP. Work on these problems is in progress. In this dynamic approach energy of 
the quasiparticle system is conserved. In contrast known attempts to use quantum discreet color variable in classical 
simulations combining dynamic differential equations for quasiparticle coordinates and momenta and discrete color 
variables results in violation of energy conservation low due to unphysical increase of kinetic energy [J| . 

II. THERMODYNAMICS OF QGP 
A. Basics of the model 

Our model is based on a resummation technique and lattice simulations for dressed quarks, antiquarks and gluons 
interacting via the color Coulomb potential. The assumptions of the model are similar to those of [lO[ : 

I. Quasiparticles masses (m) are of order or higher than the mean kinetic energy per particle. This assumption is 

based on the analysis of lattice data J19l |20|. For instance, at zero net-baryon density it amounts to m > T, 
where T is a temperature. 

II. In view of the first assumption, interparticle interaction is dominated by a color-electric Coulomb potential, see 

Eq. ([T]). Magnetic effects are neglected as subleading ones. 

III. Relying on the fact that the color representations are large, the color operators are substituted by their average 

values, i.e. by Wong's classical color vectors (8D in SU(3)) with the quadratic and cubic Casimir conditions 

IV. We consider the 3-flavor quark model. Since the masses of the 'up', 'down' and 'strange' quarks extracted from 

lattice data are still indistinguishable, we assume these masses to be equal. As for the gluon quasiparticles, we 
allow their mass to be different (heavier) from that of quarks. 

The quality of these approximations and their limitations were discussed in [lO[. Thus, this model requires the 
following quantities as a function of temperature and chemical potential as an input: 

1. the quasiparticle masses, m, and 

2. the coupling constant g^. 

All the input quantities should be deduced from the lattice data or from an appropriate model simulating these data. 



B. Color Path Integrals 

Thus, we consider a multi-component QGP consisting of N color quasiparticles representing Ng dressed gluon and 
of various flavors Ng dressed quarks and Ng antiquarks. The Hamiltonian of this system is H ^ K + U^' with the 
kinetic and color Coulomb interaction parts 



Here i = 1,. . . ,N is the quasiparticle index, N = Ng + Ng + Ng, i and j summations run over quark and gluon 
quasiparticles, 3D vectors ri are quasiparticle spatial coordinates, the Qi denote the Wong's color variable (8D-vector 
in the SU{3) group), T is the temperature, and fig is the quark chemical potential, (Qi\Qj) denote scalar product of 
color vectors. We use relativistic kinematics, as seen from Eq. ([T]). Nonrelativistic approximation for potential energy 
is used to disregard magnetic interaction and retardation in the Coulomb interaction. In fact, the quasiparticle mass 
and the coupling constant, as deduced from the lattice data, are functions of T and, in general, fj,g. 

The thermodynamic properties in the grand canonical ensemble with given temperature T, baryon (fig) and strange 
(lis) chemical potentials, and fixed volume V are fully described by the grand partition function[28| 
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where Ng = Nu + Nd + Ng and Ng = Nu + Nd + Ng are total numbers of quarks and antiquarks of all flavours, 
respectively, Ng is the number of gluon quasiparticles, p{r,Q,<T; Nu, Nd, Ng, Nu, Nd, Ng, Ng-, (5) denotes the diagonal 
matrix elements of the density matrix operator p = exp(— /3iJ) and /3 = l/T. Here a, r and Q denote the multy- 
dimensional vectors related to spin, spatial and color degrees of freedom of all quarks, antiquarks and gluons. The a 
summation and integration drdpQ run over all individual degrees of freedom of the particles, dfiQ means differential 
of the group SU{3) Haar measure. Usual choice of the strange chemical potential is /is = (nonstrange matter). 
Therefore, below we omit ps from the list of variables. 

Since the masses and the coupling constant depend on the temperature and baryon chemical potential, special care 
should be taken to preserve thermodynamic consistency of this approach. In order to preserve the thermodynamic 
consistency, thermodynamic functions such as pressure, P, entropy, S, baryon number. Nb, and internal energy, E, 
should be calculated through respective derivatives of the logarithm of the partition function This is a conventional 
way of maintaining the thermodynamical consistency in approaches of the Ginzburg-Landau type as they are applied 
in high-energy physics, e.g., in the PNJL model. 

The exact density matrix p ~ e^^^ of interacting quantum systems can be constructed using a path integral 
approach [22, l23| based on the operator identity e~^^ = e~^^^ ■ e~^^^ . . . e~^^^ , where the r.h.s. contains n + 1 
identical factors with A/3 = (3/{n + 1). which allows us to rewrite |29| the integral in Eq. ^ as follows 



J2 / dr(")d/ig(0)p(r(0),QW,a;{7V};/3) = 

P9 p, 

dpQ^^^ /dr(°)dr(i)...dr(")p(r(°),r(i\...r(");g(°);{iV};/3), (3) 



where {N} — (Nu, Nd, Ng, Nu, Nd, Ng, Ng) is introduced for briefness. The spin gives rise to the spin part of the 
density matrix (S) with exchange effects accounted for by the permutation operators Pg, Pq and Pg acting on the 
quark, antiquark and gluon degrees of freedom and the spin projections a' . The sum runs over all permutations with 
parity Kp and up-. In Eq. ([3]) 



p''' 



(r('-i),g('-i);rW,Q«;{iV};A/3) = (rC-D |e-^^^|rW)<5(Q('-i) ~ Q^)), (4) 



is the off-diagonal clement of the density matrix. Since the color charge is treated classically, we keep only diagonal 
terms ((5(Q''^^-' — Q^'-*)) in color degrees of freedom. Accordingly each quasiparticle is represented by a set of coordi- 
nates {rl , . . . jTj } ("beads") and a 8-dimensional color vector Ql in the SU{3) group. Thus, all "beads" of each 
quasiparticle are characterized by the same spin projection, flavor and color charge. Notice that masses and coupling 
constant in each p*^') are the same as those for the original quasiparticles, i.e. these are still defined by the actual 
temperature T. 

The main advantage of decomposition ([3]) is that it allows us to use perturbation theory to obtain approximation 
for density matrices p^''. which is applicable due to smallness of artificially introduced factor l/{n + 1). This means 
that in each p^'^ the ratio g^{T, ^q){Qi '\Q- ')/4:'!T\r\ ' — rj- \T{n + 1) can be always made much smaller than one , 
which allows us to use perturbation theory with respect to the potential. Each factor in ^ should be calculated with 
the accuracy of order of l/(n -I- 1)^ with 9 > 1, as in this case the error of the whole product in the limit of large n will 

be equ al to zero. In the limit (n + 1) — ;► oo p' can be approximated by a product of two-particle density matrices p\J 



[22| - |24l | . Generalizing the electrodynamic plasma results j24l | to the quark-gluon plasma case, we write approximate p 

p(r(0),rW,...r(");gW;{iV};/3) = 



per||0^ ^'^^^11^^ det||0(")'(o)|U„ det ||0(")^(o)|U, frA^W 



r(").(o) 
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In Eq. ([5]) the effective total color interaction energy 

-I »i + l -I N 

^=;7^E^ E '^u(}rt''-rr''\Arf'-rf'lQr\Qr)- (6) 

is described in terms of the off-diagonal elements of the effective potential approximated by the diagonal ones by 
means of A^<^^i(xf-'\x^;, Q„ Q,) « A/3 \^^i ix^'\x^'^ ,Q,, Q,) + <I>^^(4\ 4\ g„ Q,)] /2 ex - In (^^) . Here 
the the diagonal two-particle effective quantum Kelbg potential is 

^n^^,^S^Q^^Q.^P)- 'J:':i^^f:^ {^---''''^''+V^-S [l-erf(x|?)]}, (7) 



47r|rr - r^ 



„(0 _ i^(') ^(Oi 



„ Ir; ' — rj l/AAy, AA^- = ^J^jFK^/rriij , rriij = rriimj / {rrii + rrij) Other quantities in Eq. ^ are defined 



as follows: A^ — Ai^\/0.57r/(/3m)^ with Aa = A/27r/3/mo being a thermal wavelength of an a type quasiparticle 
(a = q,q,g). The antisymmetrization and symmetrization takes into account quantum statistics and results in 
appearing permanent for giuons and determinants for quarks/antiquarks. 

Functions (j)^^ = K2{zf^)/{z^'^'^)'i [zf^ = Al3m,{T, Hg)J I + ^f' /A^2 ^j-e defined by modified Bessel functions. 

Gluon matrix elements are (p. .' = ^2(2:1" j/i^ij V^tiQi — Q] ), while quark and antiquark matrix ele- 
ments (i>\y = K2{zlj'^ j/i^Lj j'^^aiQi ^ Qj )^fijj^cri,crj depend additionally on spin variables at and fia- 
vor index fi of the particle, which can take values "up", "down" and "strange", S„.,j. and 6f-j. are the Kro- 
necker's deltas, while 5e is the delta-like function of color vectors. These functions allow to exclude the Pauli 
blocking for particles with different spins, flavors and colors. Here arguments of modified Bessel functions are 
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^(n,,yu, ^ A/3m,(T,p,)Wl+ rf - r]"^ /A/32. The coordinates of the quasiparticle "beads" rf^ ^ rf^ +yf\ {I > 0) 

are expressed in terms of rl and vectors between neighboring beads of an i particle, defined as yl = J2k=i Q ' 
while Q , . ■ . , Cl" are vector variable of integration in Eq. dSj . 

The main contribution to the partition function comes from configurations in which "size " of the quasiparticles 
cloud of 'beads' is of order of Compton wavelength Ac = l/m{T, Hg) providing spatial quasiparticle localization. So 
this path integral representation takes into account quantum position uncertainty of quasiparticles. In limit of large 
mass spatial quasiparticle localization can be made much smaller than average interparticle distance. This allow 
analytical integration over 'beads' position by method of steepest decent. As result the partition function is reduced 
to its classical limit with point strongly interacting color quasiparticles. 

In the limit of n — >■ oo functions (pl^ describe the new relativistic measure of developed color path integrals. This 
measure is created by relativistic operator of kinetic energy K ~ y^p'^ + m?(T, fig). Let us note that in the limit 
of large particle mass the relativistic measure coincide with the Gaussian one used in Feynman and Wiener path 
integrals. 



III. SIMULATIONS OF QGP 



To test the developed approach we consider the QGP at zero baryon density with equal flavor fractions of quarks 
and antiqiarks (N.^ ~ N^ ~ Ng ~ N/3 = TY^ = iV^j = Ng = N/3). Ideally the parameters of the model should 
be deduced from the QCD lattice data. However, presently this task is still quite ambiguous. Therefore, in the 
present simulations we take only a possible set of parameters. The phenomenologic QCD estimations |25| of coupling 
constant, i.e. a^ = 5^/(47r), used in the simulations is displayed in the left panel of Fig. [T] The T-dependence of 
quasiparticle mass used in this work is also presented in Fig. [1] (left panel). The right panel presents the calculated in 
grand canonical ensemble the temperature dependences of the interparticle average distance r^ (Wigner-Seitz radius 
r^ = (3/47rn), n is the density of all quasi particles). The quark quasiparticles degeneracy parameter x and the plasma 

<T,qX^, T = 4^y J- , where 92 the quadratic Casimir value averaged over 
q2 = iv~ — 1 is a good estimate for this quantity. The plasma coupling parameter 
is a measure of ratio of the average potential to the average kinetic energy. It turns out that F is larger the unity 
which indicates that the QGP is a strongly coupled Coulomb liquid rather than a gas. In the studied temperature 



coupling parameter F are defined as: x = 
quarks, antiquarks and gluons, 02 = N^ 
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FIG. 1: Left paneliTemperature dependences of the model input quantities: the coupling constant as (scaled by 5) and 
quasiparticle mass-to-temperature ratio (m, = m,). Right panel: Temperature dependences of the calculated interparticle 
average distance r^ (Wigner-Seitz radius) , the quark quasiparticles degeneracy parameter x and the plasma coupling parameter 
r. Lines are smooth interpolation between the CPIMC points. To ilustrate sensibility of this model to input data the green 
solid symbols and arrows show changes of the CPIMC average interparticle distance (right panel) related to the quasiparticle 
mass variation on left panel. 



range, 190 MeV <T < 600 MeV, the QGP is, in fact, quantum degenerate, since the degeneracy parameter x ~ nqX^ 
(where the thermal wave length A^ was defined in the previous sect.) varies from ~ 20 to '--^ 8, see Fig. [T] (right panel). 
The degeneracy parameters for different species does not coincide, since the quasiparticle masses are different. 

As follows from analysis of Fig. [T] and Fig. [5] this model is highly sensitive to quasiparticle mass variations. This 
allow to fit lattice EOS and to choose optimal values of gluon and quark masses in order to proceed in predictions 
of other properties concerning the internal structure and in the future also non-equilibrium dynamics of the QGP. 
High sensitivity of EOS to mass variation adds to background of this procedure. Let us also note that increase 
of quasiparticles masses reduces the influence of quantum effects, so high sensitivity of EOS to masses proved the 
importance of quantum effects at these temperatures. Figure[2]presents also the trace anomaly {{e—3P)/T^) calculated 
as related combination of derivatives of equation of state. In order to avoid the numeric noise, the derivatives of a 
smooth interpolation between the CPIMC points of equation of state (Fig. [2]) were taken. The CPIMC results are 
compared to lattice data of [2|, |3| . It is not surprising that agreement with the lattice data for trace anomaly is also 
good, since it is a direct consequence of the good reproduction of the pressure. 

Let us now consider the spatial arrangement of the quasiparticles in the QGP by studying the pair distribution 
functions (PDF) gab{R)- They give the probability density to find a pair of particles of types a and 6 at a certain 
distance R = \Ri — i?2| from each other and are defined, for example, in canonical ensemble as 

gabilRi ~ R2\) ^ Z{^n\p) E/'^^^Q'^(-^i"^?)'^(-^2-r2')p(r,g,a;iV„,JVrf,iV3,iV„,iVrf,JV„JVg;/3), (8) 

The PDF depend only on the difference of coordinates because of the translational invariance of the system. In a 
non-interacting classical system, gab{R) = 1, whereas interactions and quantum statistics result in a re-distribution 
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FIG. 2: Equation of state (left panel) of the QGP from CPIMC simulations (open squares). The solid line is a smooth 
interpolation between the CPIMC points. The solid symbols illustrate the changes in equation of state for presented by Fig. [l] 
quasiparticle mass variation. Trace anomaly (right panel) of the QCP from CPIMC simulations is presented by solid line. Both 
CPIMC results are compared to lattice data of [2, 13 

of the particles. At temperatures T = 193 MeV the PDF averaged over the quasiparticle spin, colors and flavors are 
shown in the panels of Fig. |31 
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FIG. 3: Pair distribution functions of identical (left panel) and different (right panel) quasiparticles at temperature T = 
193 MeV. The distance is measured in units of cr = 1.1 fm. 

At distances, R/a > rs{rs ~ 0.35) all PDF are practically equal to unity (Fig. [3]) like in ideal gas due to the 
'Debye' screening of the color Coulomb interaction. A drastic difference in the behavior of the PDF of quarks and 
gluons (the anti-quark PDF is identical to the quark PDF) occurs at distances R/<t < r^ . Very well pronounced is the 
short-range structure of nearest neighbors. Here the gluon-gluon and gluon-quark PDF increases monotonically when 
the distance goes to zero. Behavior of these PDF is a clear manifestation of an effective pair attraction. Attraction 
suggests that the color vectors of nearest neighbor quasiparticles are anti-parallel. The QGP lowers its total energy 
by minimizing the color Coulomb interaction energy via a spontaneous "anti-ferromagnetic" ordering of color vectors. 
From physical point of view this means formation of the gluon-gluon and gluon-(anti)quark clusters, which are more 
or less uniformly distributed in space. The last conclusion comes from the fact that the (anti) quark- (anti) quarks PDF 
are close to unity. Oscillations of the PDF at very small distances of order R < O.ltr are related to Monte Carlo 
statistical error, as probability of quasiparticles being at short distances quickly decreases. 

The main physical reason for arising the well separated from each other the gluon-gluon and gluon-(anti)quark) 
quasiparticle clusters is the consequence of the great difference in values of quadratic Casimir invariants q2 ~< 
Qi \Qi > for gluon-gluon q2 = 24, (anti)quark-(anti)quark (72 = 4 and gluon-(anti)quark (72 = \/24 * 4 « 10 in color 
Kelbg (Coulomb) potentials of quasiparticles interactions . 

The second physical reason of the PDF difference is spatial quantum uncertainty and different property of Bose 
and Fermi statistics of gluon and (anti)quark quasiparticles. Fermi statistics results in effective quark-quark and 
aniquark-antiquark repulsion, while Bose one results to effective qluon-qluon attraction. More strong interaction of 



the QGP quasiparticles with gluons than with quarks is also consequence of the gluon better locaUzation in comparison 
with quarks. As we mentioned before uncertainity in particle localization is defined by the ratio l/in. Localization 
is better for heavier gluon quasipartcles. To verify the relevance of all these trends, a more refined color, flavor, and 
spin-resolved analysis of the PDF is necessary, together with simulations in a broader range of temperatures which 
are presently in progress. 

Details of our path integral Monte-Carlo simulations have been discussed elsewhere in a variety of papers and review 
articles, see, e.g. [23l . [2a | and references therein. The main idea of the simulations consists in constructing a Markov 
chain of quasiparticle configurations. Additionally to the case of electrodynamic plasmas here we randomly sample 
the color quasiparticle variables Q according to the group SU(3) Haar measure. 

IV. CONCLUSION 

Color quantum Monte Carlo simulations based on the quasiparticle model of the QGP are able to reproduce the 
lattice equation of state and other thermodynamic functions (even near the critical temperature) and also yield 
valuable insight into the internal structure of the QGP. Our results indicate that the QGP reveals quasiparticle 
clustering and liquid-like (rather than gas-like) properties even at the highest considered temperature. 

For our simulations we have introduced new relativistic path integral measure and have developed procedure of 
sampling color quasiparticle variables according to the group SU(3) Haar measure with the quadratic and cubic 
Casimir conditions. Our analysis is still simplified and incomplete. It is still confined only to the case of zero baryon 
chemical potential. The input of the model also requires refinement. Work on these problems is in progress. 

We acknowledge stimulating discussions with Prof. P. Levai, D. Blaschke, R. Bock, and H. Stoecker. 
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